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Abstract 

Let f(x,p) be a function on the phase space. The function / 
corresponds to the classical observation. Let p(x,p,£) = <p 2 {x,p,£) be 
some non-negative density function on an extended phase space. The 
density function p corresponds to the generalized state on the extended 
phase space. Then, as usual, the observation / in the generalized state 
p is the value < /, p >=< f, ip 2 > of the integral of the function / in 
the distribution p. It is supposed, that in quantum observations all 
distributions p are realized not, and only distributions p = ip 2 , where (p 
belongs to some linear subspace H. (averaging wave functions) in space 
of all functions on the extended phase space. Besides it is supposed, 
that in quantum experiments values of the spectrum of the linear 
operator Aj of the quadratic form < /, tp 2 >=< ip,Afip 2 >, where 
(p £ 7i, are observed. In this paper we consider certain hypotheses on 
the averaging process of wave functions. Then it is shown that the 
spectrum of usual operator of the quantum observable corresponded 
to / is small differing from the spectrum of the operator At. Other 
side, this approach enables one to recover probability distributions in 
the phase space for wave functions. 



1 



Introduction 



In this paper there is considered a mathematical model of microworld pro- 
cesses based on certain constructions of the geometric quantization [Q which 
takes in account the fluctuating affects of the media, such as fluctuations of 
the vacuum. We show that an enlargement of the phase space and of its mo- 
tion group and an account for the diffusion motions of microsystems in the 
enlarged space, the motions which act by small random translations along 
the enlarged group, lead to observable quantum effects. 

I show that certain conventional nonrelativistic descriptions of quantum 
systems can be obtained as asymptotic approximations to the model sug- 
gested in this paper with respect to the power series expansion in Planck's 
constant. The parameters of the proposed model are estimated on the base 
of the experimental data of Lamb's shifts in the spectrum of hydrogen's 
atom [|LR|1 . In quantum mechanics these shifts can not be explained in the 
framework of the nonrelativistic quantum model and the idea of vacuum's 
polarization around a point source (radiational corrections) is incorporated 
to do the trick. 

On the one hand, the suggested model can be used to deepen our per- 
ception on the nature of quantum effects; on the other hand, it leads to yet 
another construction of quantizations of mechanical systems, the construc- 
tion that might be of help in some exceptional cases. 

Here I propose to consider probability amplitudes on the enlarged phase 
space which includes the space of inner states of the system, its own in each 
experiment. I assume that in the experiments on quantum systems our clas- 
sical gadgets can only register the probability distribution corresponding to 
the averaged probability amplitudes over fluctuating actions on the quantum 
system. I will show that in the suggested model the averaged probability 
amplitudes on the phase space are parametrized by complex functions on the 
configuration space, i.e. the wave functions. In particular, formula (|8]) of 
Teorem 4.1 obtained in this model reflects the probability distribution for a 
wave function in the phase space. Wigner ||Wfl who obtained so-called "qua- 
sidistributions" in the phase space was first solving this problem. However in 
certain cases his quasidistributions can be negative (unlike the ones obtained 
here) and therefore have no physical meaning. 

In sec. 1 the main hypotheses on quantum observations are formulated 
and a formulation of the main problem is given. 
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In sec. 2 the action of the enlarged group of motions on the enlarged 
phase space is described. 

In sec. 3 there are introduced hypotheses on the averaging operation and 
the theorem on the form of averaged distributions in the enlarged phase space 
is formulated. 

In sec. 4 the main result of the paper is proved: the asymptotic proximity 
of linear operators of the considered model to the linear operators of quantum 
observables and the probability distribution in the phase space for a wave 
function. 

In sec. 5 there are listed several unsolved problems and trends of further 
investigations. 

Appendix contains proof of the theorem on the form of averaged distribu- 
tions in the enlarged phase space and an estimate of the model's parameters 
based on the experimental data of Lamb's shift in the spectrum of hydrogen's 
atom. 



1 Main hypotheses and formulation of the 
problem 

First let us recall the main notions of classical mechanics for the flat config- 
uration space, cf. e.g. [DFN, |B2"|. Let M (locally, R 2n ) be the phase space 



with coordinates of its points being (q,p), where q = (g 1; q n ) is the point's 
position in the configuration space and p = (pi, ...,p n ) is its momentum. In 
classical mechanics the dynamics is given by the formula df/dt = {f,H}, 
where {., .} is the Poisson bracket determined by the exterior 2-form u (lo- 
cally of the form J2i<i< n dqiAdpi) and if is a fixed function, the Hamiltonian. 

The canonical transformations g of the phase space M are (smooth, 1-1) 
maps g : M — > M preserving to. Let S denote the group of all canonical 
transformations (symplectic diffeomorphisms) of the phase space. The dy- 
namics of a mechanical system is given by a 1-parameter subgroup of S. 
To observables of a mechanical system the functions / on the phase space 
are assigned and to the states there are assigned nonnegative distributions 
pdqdp, or shortly, p. The mean value of an observable / in the state p is by 
definition 

(f,p) = / f(q,p)p(q,p)dqdp. 
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We will make the following hypotheses on a quantum system on the above 
phase space: 

1. The phase space is enlarged to the space E which is the total space of a 
fibration 7r : E — > M over M with fiber F which is a manifold of 'inner 
states' of the system (if M = R 2n then E = R 2n x F). 

2. Let pr : P — > S be the nontrivial central extension of the group S with 
the help of circle T = R/hZ (Lie algebra of P is the Poisson algebra). 
The group P acts by diffeomorphisms on E so that n o g' = pr(g') o n 
for any g' e P where g' : E — > E. 

3. The observables of a quantum system, as well as of a classical one, are 
given by functions on the phase space. 

Denote by p(q,p, ^)dqdpd^ a nonnegative distribution on E determined by 
a distribution density p(q,p, £) with £ e F and a P-invariant measure dqdpd^ 
on E. Let us present p in the form p(q,p,£) = (p 2 (q,p,£)i where (p(q,p,£) is 
a function on E such that 



Denote by p(q,p,£) = <p 2 (q,P,0 an averaged density distribution and 
<p(q,p, obtained from ip(q, p, £) by an averaging process under the action of 
small fluctuations of the system on E (a mathematical model of this averaging 
process will be given in sec. 3. 

4'. Only averaged density distributions of the form p(q,p,£) = (p 2 (q,p, £) 
are realized in quantum measurements. The map <p i— > <p is a linear projection 
operator whose form is given in sec. 3, where the refined hypothesis 4), is 
given. 

5. Only averaged values of observables f of the form 



are realised in measurements, where p(q,p,£) = <p 2 (q,p, £) and (p is the av- 
eraging function on E. 

Now, consider the completed Hilbert space H C L 2 (E , dqdpd^) of aver- 
aged functions <p(q,p,£) on E with respect to the standard inner product 
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To a function / on the phase space assign the linear operator Af in 7i 
given by the formula 

&,A^) d = f (f,p) = (f,^). 

The main result of this paper is the proof of the fact that under the 
natural P-action on E the operators Af are approximately (up to terms of 
order h, where h is Planck's constant) coincide with the operators of quantum 
observables in the accept definition of quantum mechanics. 

2 A description of the group P and its action 
on the enlarged space of states 

Let pr : P — ► S be the nontrivial central extension of the group S of 
canonical transformations with the help of the circle (1-dimensional torus) 
T = R/hZ. Let po(2n) be the Poisson algebra, the nontrivial central exten- 
sion of the Lie algebra h(2n) of Hamiltonian vector fields. In a sense that 
Lie theory is applicable to infinite dimensional case po(2n) and h(2n) are 
the Lie algebras of the groups P and S, respectively. 

Two bundles E, E' over a symplectic manifold M (with fibers F, F', re- 
spectively) with P-actions on them are equivalent if there is a diffeomorphism 
of bundles with P-action compatible with projections and the S-action on 
the base. To describe bundles with P-action over M we have to consider first 
the structure of the fiber F of the bundle E in detail. 

Denote by So the subgroup of the canonical transformations S that pre- 
serve the origin, let Po = pr _1 (So). Clearly, P = T © So. 

Since P C P acts on E and preserves the fiber F over the origin, a 
P-action in F is defined. In particular, T C P acts in F. 

Theorem 2.1. ([K]). A bundle E = R 2n x F with a P-action compatible 
with an h(2n) -action on the space of functions on R 2n is uniquely up to 
equivalence of P -bundles determined by the space F with a P -action. 

To a P-action on E we assign in the standard way a Lie algebra homomor- 
phism po(2n) — > vect(P), the derivative of the P-action on E. Denote by 
t the vector field on F corresponding to the action of the 1-parameter group 
generated by T on E (r can also be defined as the value of the derivative of 
the action T — T t with respect to the parameter t at t — 0). 
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Corollary 2.2. In any P -bundle E one can chose a trivialization (iso- 
morphic to R 2n x F) so that the vector fields on E corresponding to Hamilto- 
nians linear in q and p, i. e. of the form H = Y,?=i( x iPi ~ UiQi) + c G po(2n), 
are given in local coordinates (q,p,£) on E by the expression 

n n 

Dh = ^{xidqi + yidpi) + (c - ^ yiqi)T. 
i=i i=i 

Equivalently, in the global form, this means that for an arbitrary function 
<f(q,p,C) 071 E the action of the one-parametric subgroups Gf for the above 
H G po(2n) is given by the formula 

G?(p(q,P,Z) = v{q + tx,p + ty,T t(c ^ {yyq)) {£)), 

where (y, q) = Y%=i UiQi an d Tt( c -{y,q)) ^ the element of the 1-parameter group 
T corresponding to the value t(c — (y, q)) of the parameter. 
In what follows we set: 

W t x ' v = G? (or H = j2(x lPl -y iqi ). 

i=l 

Denote by W the subgroup of P generated by Wf ,v for (x, y) G R 2n (usually 
W is called Heisenberg-Weyl group). Clearly, W is the inverse image of R 2n 
with respect to the projection pr : P — > S, where R 2n C S is considered as 
the subgroup of translations. 

3 Averaging 

Now, let us pass to averaging of a function <p(q,p,£) on a P-bundle E. We 
need it in hypothesis 4. 

Let us start with the assumption that the averaging if \— > <p is associated 
with a diffusion process (Brownian motion) on a P-bundle caused by inac- 
curacy of the setting of the quantum observational device. This process acts 
locally by small shifts by Aq and Ap along the coordinates q and p of R 2n , 
respectively, and globally with the help of W^ q ' Ap G W acting on E. 

More precisely, let r be the diffusion time. The functions are transformed 
as for a diffusion process, cf. e.g. [I]: 

V9(r + Ar,g,p,0 = J K(Aq, Ap, Ar) x 

R 2n 
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xW^Mr,q,P,Od(Aq)d(Ap) + o(Ar), (1) 



where K(Aq, Ap, At) is the probability density of shifts by the vector 
(Aq, Ap) during the time At. 

Naturally, we make the usual assumptions about K(Aq, Ap, At): 

— K(Aq, Ap, At) is rapidly decreasing at infinity; 

- the mathematical expectation of the shift vector (Aq, Ap) is zero; 

- the diagonal elements of the matrix of 2nd moments are of the form 



where a and b characterize the "intensity" of shifts along positions and mo- 
mentum; 

- shifts along distinct directions are poorly correlated with each other, 
i.e. the offdiagonal 2-nd moments are of order o(At); 

- the moments of K(Aq, Ap, At) of orders greater than 2 are also of 
order o(At). 

Starting from ([!]) and the above assumptions we expand ip{r,q,p,^) in 
the Taylor series to derive as At ^ a differential equation similar to the 
diffusion equation. The equation represents the averaging of the function (p 
over fluctuations. The asymptotic of solutions of this equation as r — > oo 
is given by the following theorem, where j is the imaginary unit, * denotes 
the complex conjugation and ip(x,^) is a complex- valued function on R n x F 
such that 



Theorem 3.1. Lei <p(r, q,p, f) satisfy ffi) and ip(0, q,p, f) = <p(q,p,£), 
where f£ ip(q,p,T t (£))dt = 0. Then (p(r,q,p,£) asymptotically approximates, 
as t — > oo , to the function 




2a, 2 Ar + o(Ar), 




2b\ At + o(At) 






(3) 



where 




4 / b x ...b. 



a\...a. 



n 



■n 



) 
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x L(x, exp (-j 2 -^~\ + r (x, exp ) dx. (4) 

The function ij)(x,£) is obtained from the function (p(q,p,£) according to the 
formula 

( it " hi, . 2 \ ( .2n(p,x)\ , , , 
xexp^-- 2^ — (ft - Xi) Jexp^j j dqdpdt. (5) 

Besides, iftp{x,^) is an arbitrary complex-valued function satisfying (^) then 
the composition of transformations if) i— > (p and (p t— > ip given by ([Jj and (^) 
is identity. 

Proof is given in Appendix. 

Now we are ready to refine hypothesis 4' from sec. 2 as follows: 
4. Let (f(q,p,£) be a real-valued function on E such that 

[ h <p(q,p,T t (O)dt = 
Jo 

and p(q,p,£) = f 2 {q,P,0 be the probability density on E. The averaging 
operation mentioned in hypothesis 4' from section 2 is caused by the fluctu- 
ation process described by equation ([Jj where ip is the initial state and <p is 
the asymptotic one as the time of fluctuation tends to infinity. 

A corollary of this hypothesis: the map ip i— > (p is obtained (by Theo- 
rem 3.1) as the composition of maps given by formulas ([5]) and (f|) respec- 
tively. 

Consider the Hilbert space TC of real functions (p on E of the form (Q) 
and with the standard inner product 

Up>, (p") = f &fl'dqdpd£. 
Je 

Denote by Ti the Hilbert space of complex- valued functions ib(x,£) on R n x F 
satisfying (|2|) and with the inner product given by the formula 

tyxM = Re il>i(x,£)if>Z(x,£)dxd€. 

JR n xF 
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Making use of Theorem 3.1 we directly derive the following 

Corollary 3.2. The map given by formula (fjj is an isomorphism 
of Hilbert spaces H and H. 

A function ip(x,^) satisfying (|2|) will be called a wave function, the cor- 
responding function <fi{q,p,!;) will be called the probability amplitude in the 
enlarged phase space. 

Remark. Formula (|5]) gives the so-called position representation of 
ip. The momentum representation of ip is given by the Fourier transform of 
?/>(x, £) with respect to x. We will not ned it. 

4 An estimate of the operator of an observ- 
able 

In this section we will estimate the operator Af acting in the Hilbert space 
"H. of functions of the form (^). By definition Af is given by the following 
expression: 

(f,p)n= (0, A f0)-H= f(q,p)fi 2 (q,p,O d( i d P d C ( 6 ) 

J E 

Since by Corollary 3.2 the Hilbert spaces 7i and 7i are isomorphic, we 
can estimate Af by estimating image of Af uiTi under this isomorphism. 
Let us substitute (f|) in (|j). Represent the function 

V 2 (q,P,0 = (J B(q,p,£,x)dx^J , 
where B(q,p,^,x) is the integrand of (^), in the form 

ft 2 (q,P,0= / B(q,p,£,x')B(q,p,£,x)dx'dx. 

JR n JR n 

Since the integrals of the product ip(x', C)ip{x, £) and of the product 
ip*(x', £)ip*(x, £) over £ are zero due to (0), we get after simplification 
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X 



/// /(9>P) ex P 



R n R" E 



n i=i a 



x^(x',^)^*(x,C)exp -j 



27r(p, x — x') 



dqdpd^dx'dx. 



This implies that the kernel of Af on the space 7i is of the form 
2\i fh...b n Y 



A f (x,x') 



a±...a n 



h 3 

x / /(g.p)exp 

27r(p, x — x') 



i? 2 

x exp -j 



7T " 6 
i=l 



dqdp, 



(7) 



and the density of the probability distribution p(q,p) = v5 2 (<?,p) correspond- 
ing to i[j(x,£) is given by the next theorem. 

Theorem 4.1. Let ip(x,£) be a wave function (in the position 
representation) on R n x F satisfying (fjj; then the density of the density of 
the probability distribution on the phase space is given by the formula: 



x 



_2_\i / h...b n 
h 3 J \ai...a n 

exp 



F R n R n 

x exp —j 



h i=l a 

2n(p, x — x') 
h 



ifj(x , ,^*(x,^)dx'dxd^ (8) 



This distributions is different from Wigner's quasidistributions by the 
integration with the function 



2\ t (b x ...b T , 
h 



d\ . ..Qj r 



exp 



7T 



n b 



7E>-^) 2 + fe-^) 2 ) 

n i=l a 



Now in order to represent the operator Af in the form habitual in quan- 
tum mechanics (cf. [ KYa ]), substitute in ([7]) the expression of / in terms of 
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its Fourier transform: 

f{q,p) = (^-) J f{u,v)exp[-j({q,v) + (p,u))]dudv, 

where n 

/(g,P) = (^j J f{u,v)exp\j((q,v) + (p,u))}dqdp. 

R 2n 

In the obtained formula we first integrate over q, then over p and u. We 

get: 

1 r ~ (2tt{x-x') \ ( . Ix + x' \\ , 

= ^i h { — —> v )^\- 3 \—> v )r> 

where 

f h {u,v) = /(«,«) exp [~t{^ + V)) • 

If instead of fh{u,v) we take its Taylor series expansion in powers of h 
up to order n we get the corresponding asymptotic representation of Af. In 
particular, the asymptotics of the 0-th term of the expansion is 

A ' {X ' X "> = hU ») exp [-3 «-)) *. 

This expression coincides with the expression for the operator of an ob- 
servable in quantum mechanics (formula (14) in [FYa] whose h is our h/2iv). 
This implies the next result of the paper: 

Theorem 4.2. The linear operator Af of the form (Q) in the Hilbert 
space 7i which is constructed from the classical observable f under assump- 
tions 1-5 on the process of quantum observation is asymptotically close to the 
conventional operator of quantum observable given in the position represen- 
tation. 
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If a more precise estimate of Af is required we can always take into 
account more terms of the Taylor series expansion of fh(u, v) with respect to 
h. 

Examples of exact formulas for Af. By the usual abuse of language let 
us denote the operator of multiplication by a function F by F. We have: 

Aq i Xj, 

. h d 



2w dxi 
A q * = x 2 + / '"' 



h 2 d 2 hh 

A n 2 = — — 77 + 



P * 47T 2 dx 2 4:710,1 ' 

It follows that for the Hamiltonian / = p 2 / (2m) + muj 2 q 2 /2 of the linear 
oscillator with eigen frequency u we have 

h 2 d 2 muj 2 x 2 Jt, h(b 2 + m 2 uj 2 a 2 ) 



8ir 2 m dx 2 2 r-J 8iraibim 

which differs from the conventional operator of the quantum linear oscillator 
by constants (the last summand). 

In other side, if / = f(q), i.e. does not depend on momenta, we deduce 
from (0) after integration over p and x' that Af is the operator of multipli- 
cation by 



The passage / i— > / is the convolution of / with the density of the prob- 
ability distribution with dispersion along the g^-axis equal to haiK^bi). 
Suppose 

ai/bi = a 2 /b 2 = a 3 /b 3 = a/b. 
Then for the Hamiltonian 

f(q,p) = (Pi+P2+Pl)/( 2m ) + V(qi,q 2 ,qa), 
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we have 

h 2 ( d 2 d 2 d 2 \ 3hb -. 

In particular for the Hamiltonian of the hydrogen atom whose 

V{qi,q 2 ,qz) = -e 2 /r, (10) 

where e is the charge of the electron and r 2 = J^^ilj, the operator Af 
differs from the operator by an irrelevant constant 3/i6/(47ra) and the extra 
smoothness of the Coulomb potential. Thus our hypotheses predict that 
the theoretical spectrum of the hydrogen atom computed on the base of the 
conventional Hamiltonian should differ from that obtained in experiments. 
Such a descrepancy of theory and experiment was indeed detected in 40s 
[LR]. It is called Lamb's shift of hydrogen atom's levels and is conventionally 
explained in quantum electrodynamics by an interaction of the election with 
a fluctuating electromagnetic field. Comparison of these experimental data 
with calculations of the spectrum of Af given by (§), (|T0"1) via perturbation 
theory yields the following estimate of parameters of our model (details see 
in Application A. 2): 



a/b = 3,41- 10 4 sec/gr; Aq = y — = 4, 24 • 10~ 12 cm. 

Hence the standard deviation of the normal distribution Aq is comparable 
with h/(27cmc) = 3, 8 • 10 _11 cm, the minimal position error of electron (in 
the rest frame) obtained in quantum electrodynamics. 



5 Conclusion and plan of further study 

The description of quantum systems is usually based upon certain formal 
procedures starting from a classical description of the corresponding mechan- 
ical systems. Many physicists and mathematicians, starting with Einstein, 
searched for a meaning of these procedures but the success of quantum me- 
chanics approved the formal approach to the quantization procedure and von 
Neumann's theorem on "hidden parameters" [Nu] discredited for a long time 
such a search as a direction of scientific investigations. 
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The interest to the problem was revived in 50s in works of Bohm and de 
Broglie [B], [Br] and maintained in a number of later works, cf. [BV], [N], 
[PG], [M], [Ba], [KV]. At the same time the rigidity and lack of motivation in 
the mathematical requirements in von Neumann's theorem became manifest 

A detailed analysis of the problem of introduction of hidden parameters in 
quantum mechanics and von Neumann's theorem is given in [Kh]. In partic- 
ular, there is given a formal model with hidden parameters for a "solitary" 
quantum system. The difficulty of introducing the classical probabilistic 
model for quantum phenomena is associated with a non-local character of 
these phenomena, confirmed in a number of experiments, cf . the review [Gr] , 
[SM] (however, cf. [Bl], Appendix 1). 

We have shown how to deduce operators of quantum observables on the 
base of hypotheses of sec. 1. The reason causing fluctuations on the ex- 
tended phase space are diverse: a fluctuating external force, inaccuracy in 
the description of the real system, etc. We have shown that irrespective of 
the nature of the fluctuation quantum effects will be observed in such sys- 
tems with accuracy determined by the Planck's constant and the ratios Oj/fei, 
where dj and hi are the intensities of random shifts along the i-th position 
and momentum, respectively. 

Elsewhere I intend to investigate the following problems: 

- Take into the account relativistic effects. We have constructed Af 
by formula (0) having given a classical observable / under the assumptions 
that the intensities of shifts a and b along the positions and momenta are 
constants. This is not a relativistic requirement and it is desirable to find the 
dependence of the intensities on momenta and the masses of particles and to 
refine formula (§). 

- Take into the account spin. For this we should replace our group P by 
a supergroup whose Lie superalgebra is po(2n|m) or its "odd" counterpart 
b A (n), see [L]. 

- Generalize our construction to phase spaces more general than direct 
products of the configuration space by the space of momenta. This should 
lead to restrictions on the Planck's constant h cases by the geometry of the 
phase space, cf. [B2]. 

- Describe dynamics of observable taking into account fluctuating ac- 
tion of the ambient media. As is shown in App.l the time for stabilization 
(averageing) over fluctuations is of order h/(2irab). If this quantity is small 
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the dynamics can be approximately described as the superposition of a fast 
and a slow movements. The fast motion leads to averageing of probability 
amplitudes and the slow one describes the motion in the space of averaged 
amplitudes. The classical Schroedinger equation only describes the slow con- 
stituent of the motion of a microobject. 



Since (3.1) is determined via the action of the Heisenberg-Weyl group Wf' y in 
the space of square integrable (and complexified for convenience) functions on 
E, let us decompose a function (p(r,q,p,£) into the integral over irreducible 
representations of W. 

The decomposition will be performed into three steps. First, let us con- 
sider the action of T = R/hZ in the space of functions on F. Let </?(£) be a 
function on F. Since ip(T t (x)), where T t G T, is periodic in t with period h, 
it has the Fourier series expansion 



In particular, for t = we get </?(£) = J2kez <Pk(Q- It follows from the 
definition of <Pk(Q that </?(£) i— > v?fc(£) is a projection and T acts on the 
functions v?fc(£) by the formula 



As is easy to see, the functions tpk and tp n are orthogonal for n ^ k. 
A function ip(q, p, ^) on the P-bundle E = R 2n x F can also be represented 
in the form 



Appendix 1. Proof of Theorem 3.1 




where 





(Al) 



(A2) 



where 
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Having substituted (A. 2) into (3.1) and taking (A.l) into account we get the 
following equations for each of the orthogonal components of tp in coordinates 
on E given by Corollary 2.2 

ip k (T + AT,q,p,0= f K(Aq,Ap,AT)(p k (T,q + Aq,p + Ap,Z) 

x exp ^_j 2irk ( A P> g) ^ d(Aq)d(Ap) + o(Ar). (A3) 

For k — the equation (A.3) turns as At — > (as is clear from the Taylor 
series expansion) into the following differential equation on (po(r, q + Aq,p + 
AM): 

C^q_ " / 2 gVo , 2 gVo\ 

^ -fer ^ + ^ ^ j' 

and since by hypothesis of Theorem 3.1 



(0,q,p,£) = 1/h j\{q,p,T t (£))db = 0, 



then (p (T,q,p,^) = for any r > 0. 

For k 7^ let us expand ^(r, q,P,£) into the Fourier integral with respect 
to p, i.e. represent it in the form: 

n 

¥k{r,q,p,0 = (^jpj J Rn 0k{r,q,x,£)exp ^j 27lk ^' x ^ d x , (AA) 
where 

n 

<p k (r, q, x, = <p k (r, q, p, g) exp ^-j 2nk ^ x ) ^ ^p. (A5) 

Remark. As is easy to verify, the functions of the form 

y3fc(r, g, x, £) exp (j2irk(p, x)/h) for fixed A;, £ (A6) 

span a ^-invariant subspace in the space of functions on E with the In- 
action described in the Corollary 2.2. 

Since the spaces of functions (A. 6) are orthogonal to each other for dis- 
tinct k,x and are ^-invariant, then equations (A.3) split into a system of 
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equations, each for each subspace, each obtained from (A. 3) by substituting 
ip k (r, q, x, f ) exp (j2nk(p, x)/h). 

As a result of all these transformations we get 



(fk(r + Ar,q,x,£) = / K(Aq, Ap, AT)(p k (r,q + Aq,x,£) 

JR 2n 



X 



x exp ^k(Ap^,x g) ^ rf(Ag ) rf(Ap ) + ( Ar)> 

which by expansion of the function 

<£ fc (T, g + Aq, x, f ) exp (j2irk(Ap, x-q)/h) 

into the Taylor series in powers of Ap and Aq and with the above properties 
of K{Aq, Ap, At) take as At — > the form 

where y3fc are functions in r, g, x, £ for G Z \ {0}. 

To study eqs. (A. 7), consider the corresponding eigenvalue problems: 

This is a stationary Schroedinger equation for harmonic oscillations [FYa]. 
If (fik(q) —> as \q\ — > oo, the equation has a discrete spectrum. Its eigen 
values are of the form 

W., fc „ = -E^T^(2^ + 1), (A8) 
i=i rt 

where ki,...,k n are nonnegative integers. The corresponding eigen func- 
tions <f>k,ki,...,k n (Q, x ) are products of the Chebyshev-Hermit polynomials in 
(bi/ai)^2~n]k\/h(qi - x t ) by exp {-n\k\/hY%=i(bi/a i )(q i - ^i) 2 ) • 

Since eigen functions of the Schroedinger equation are orthogonal to each 
other and constitute a complete system in the space of square integrable 
functions in q, then ip k can be represented as the series: 

J oo 

<Pk{r,q,x,£) = -= Ck,k u -,kn(T,x,Z) (Pk,ki,...,kn(q,x), (A.9) 

v z k u ...,k n =0 
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where 

Cfc,fc 1 ,...,fc n (r,x,0 = V2 J <p k (T,q,x,Z) <PkM,-,k n (q, x)dq. 

R" 

The multiple l/y/2 is taken for convenience. In particular, the maxi- 
mal eigen value A±i,o,-,o = — Yh=i ^cah/h is attained on the normed eigen 
function 



where 



r-(2\-(b 1 ...b n Y f , ( 

c±i,o,...,o = V 2 - 



ft,/ \ai. 



— J / <£±i exp [ -^X^— (* -^) 2 ] dg. 

■ a nJ I \ hf^di J 



Now, let us return to the nonstationary equations (A. 7). Having substi- 
tuted in them the expressions for <^fe,fc 1 ,...,fc„(?, x) are eigen functions of the 
right hand side of the equation, we get equations 

OO Q OO 

E u ^k,k\,...,k n - \ ^ -y 

k\,...)k n = Q k\ , . . . ,k n = 

which due to orthogonality of the functions <f>k,ki,...,k n (q, x) in the space of 
square integrable over q functions split into the system of equations 

, . . . ,k n 



dr 



for k E Z \ {0}, k\ > 0, > 0. The solutions of these equations are of 
the form 

Cfc,fc lr ..,jfc n (T, x, = c fcifclv .. )fcn (0, x, exp(A fc)fcli ... )fcn r). (All) 

Since \k,k lt ...,k n < 0, then Ck,k lt ...,k n decrease exponentially as r — > oo and 
the largest contribution to 0k(T,q,x,£) is given for not too small r by the 
terms of the series (A. 9) with the largest eigen value. 

Thus, with (A. 9), (A. 11), (A. 8) we get the following asymptotic in r — > 

oo: 

(p k (T, q, x, £) ~ -^c M ,...,o(0, x, £)<Pk,o,...,o(q, x) exp ^-r ^ j . 
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With (A. 4) we deduce that in (A. 2) every summand ipk with k 7^ 
exponentially decrease with the growth of r. Hence the largest contribution 
to if is given after a while by the terms with the largest exponent, i.e. for 
k — ±1. In other words, we have, asymptotically, 

(p ~ y?_i + (p 1 . 

Having substituted here, consecutively, (A. 4) for k = ±1, the asymptotic 
expressions for <p± 1 obtained above and expressions (A. 10) for (p±i t o,. ..,0 we 
finally get 

. 1 /2\? /6i...6„V / "27^1^ 

x /- exp Hs^ (g ^" i)2 ) x 

x ^c_i iO) ... lO (0, x, exp ^_j 2irk & x ) ^j + 

+ Cl ,o,...,o(0, x, exp (j 27rk ^ x ^ ) cfa, (A. 12) 

where c±i ; o,...,o(0, x, £) is the result of operations 

<p(q,p,t) = v(0,g,p,0 >-»• y±i(g,p,0 >-»• ^±i(?,p,0 ^ c ± i A ..., (0, x, £) 

according to the formulas (A. 2), (A. 5), (A. 10). 

It is not difficult to verify that if if is a real function 
then ci 5 o,...,o(0, x, £) = cl 10 ... i0 (0, x,0- Therefore, setting 

^(x,0 = c_i i0 ,...,o(0,x,0 

and having substituted this into (A. 12) we get Theorem 3.1. 
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Appendix 2. Estimate of the parameter a/6 of the 

considered model 



In this estimate we follow the method of [We] to justify the Lamb's shift. 
Consider the operator Af given by (0) for the Hamiltonian f(q,p) as a pertur- 
bation of the operator H for the hydrogen atom of the quantum mechanics. 
The perturbation theory implies that in the first order the increment SE n of 
the eigenvalue E n of H is of the form 

5E n = f p n (x)(V(x)-V(x))dx, 

where p n (x) = \ip n (x)\ 2 and ip n (x) is the eigenfunction of with eigenvalue 
E n and where V(x) = — e 2 / \jx\ + x\ + x\ = —e 2 /r. By definition the func- 
tion V(x) is the mathematical expectation V(x) = M q V{x + q), where q is 
normally distributed with density 

@) %xp (-|r ( * 2+ * 2+ « l) )- 

After simplification we get: 

5E n = / (M q (p n (x + q)) - p n (x))V(x)dx. 

J R 3 

Since p n (x) is smooth (infinitely differentiable), then if the standard de- 
viation of q is essentially smaller than the atom's radius we can make use of 
the Taylor series expansion of p n {x + q) to compute M q (p n (x + q)). We have 

Pn(x + q) « (1 + (q, V) + -(<?, V) 2 Pn {x). 

Since M q (qi) = 0, M (? (g i g i ) = for i ^ j and M q (qf) = ha/(±irb), then we 
get the approximate equality 

M q (p n (x + q)) = pn(x) + ^-Ap n (x), 

where A = d 2 jdx\ + d 2 jdx\ + d 2 jdx\. 
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Having substituted this M q (p n (x + q)) into the expression for 5E n and 
since A is self adjoint, we get 



ah f * dJVjx) J 

bE n = — / p n (x) }^ 2 dx. 
87r6 Jr3 r ±^ dx% 



Since for V(x) = —e 2 /r we have AV(x) = 4Tre 2 5 (x), where 5 (x) is Dirac 
delta-function, 

rr ^ ahe 2 , . 

5 n = ~w pni0) - 

For the hydrogen atom 

1 / TflG^ \ ^ 

P„(0) = |^(0)| 2 = ^( T ) . 

see e.g. [STZh] p. 342, where Ti = h/(2n), and therefore 

ahue 2 1 / me 2 \ 3 a m 3 e 8 a m 3 a 4 c 4 

where a = e 2 /(he) = 1/137 and c is the speed of light in vacuum. It follows 



6 m 3 o; 4 c 4 

In experiments of Lamb and Retherford [LR]it had been established that 
for the hydrogen atom SE 2 = 1058gHz = 1058 • 10 6 /i erg. Comparing this 
expression with our value of SE 2 , we directly get the estimate for a/b: 

a/b = 3,41 • I0 4 sec/gr. 

Accordingly, the standard deviation of q in each coordinate is 

A qi = ^ah/2b = 4, 24 • 10~ 12 cm, 
which is essentially smaller than the radius of the hydrogen atom. 
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